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Abstract: We investigate the cutoff effects in 2-d lattice 0(iV) models for a variety of 
lattice actions, and we identify a class of very simple actions for which the lattice artifacts 
are extremely small. One action agrees with the standard action, except that it constrains 
neighboring spins to a maximal relative angle 5. We fix 5 by demanding that a particular 
value of the step scaling function agrees with its continuum result already on a rather 
coarse lattice. Remarkably, the cutoff effects of the entire step scaling function are then 
reduced to the per mille level. This also applies to the ^-vacuum effects of the step scaling 
function in the 2-d 0(3) model. The cutoff effects of other physical observables including 
the renormalized coupling and the mass in the isotensor channel are also reduced dras- 
tically. Another choice, the mixed action, which combines the standard quadratic with an 
appropriately tuned large quartic term, also has extremely small cutoff effects. The size of 
cutoff effects is also investigated analytically in 1-d and at TV" = 00 in 2-d. 
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1 Introduction 



2-d O(N) models share many features with 4-d non-Abelian gauge theories. They are 
asymptotically free, have a non-perturbatively generated mass gap, and, for N = 3, even 
instantons and 6- vacuum effects. 2-d O(iV) models are integrable and have an analytically 
known exact S-matrix [1-3]. Based on this result, the exact mass gap has been extracted 
analytically [4]. This has even been extended to the mass gap m(L) in a finite periodic 
volume of size L [5-8] , which then provides exact information on the step scaling function 
introduced in [9]. Furthermore, 0(N) models can be simulated very efficiently with the 
Wolff cluster algorithm [10, 11]. For these reasons, 2-d 0(N) models are ideally suited as toy 
models for QCD, on which non-perturbative lattice methods can be tested with exquisite 
precision. Systematically controlling ultra-violet cutoff effects due to a finite lattice spacing 
a is a major objective of any lattice calculation. Symanzik's improvement program provides 
a reliable effective field theory basis for achieving this goal [12-15]. Interestingly, for 
many years the observed cutoff effects of the step scaling function in the 2-d 0(3) model, 
which seemed to be of order 0(a), were in apparent contradiction with the 0(a 2 ) behavior 
predicted by Symanzik's theory [16]. Only recently, a careful higher-order investigation of 
Symanzik's effective theory resolved this puzzle by identifying large logarithmic corrections 
to the 0(a 2 ) effects, which mimic 0(a) behavior [18]. 

While Symanzik's improvement program aims at reducing cutoff effects in a system- 
atic manner, order by order in the lattice spacing, the perfect action approach aims at 
completely eliminating cutoff effects at least at the classical level [19]. The fixed point 
action corresponding to a given renormalization group blocking transformation is indeed a 
classically perfect action, which is completely free of lattice spacing artifacts, even at ar- 
bitrarily coarse lattices. Remarkably, a practical parametrization of the classically perfect 
fixed point action, which includes a large number of terms beyond the standard nearest- 
neighbor coupling, was found to drastically reduce cutoff effects even at the quantum level 
[19]. Recently, the study of cutoff effects in the 2-d 0(3) model has been driven to an- 
other extreme by studying topological actions [20]. Topological lattice actions are invariant 
against small continuous deformations of the lattice fields. The simplest topological action 
constrains the relative angle of neighboring 0(N) spins to a maximally allowed angle 6, 
and assigns zero action to all allowed configurations. This action does not have the correct 
classical continuum limit, it cannot be studied in perturbation theory, and it violates the 
Schwarz inequality between action and topological charge for N = 3. Hence, one may 
consider this action as tree-level impaired (rather than e.g. 1-loop Symanzik improved). 
Despite these deficiencies, the 2-d 0(3) model with a topological lattice action was still 
found to have the correct quantum continuum limit [20]. Its cutoff effects at practically 
accessible correlation lengths were found to be even smaller than those of the standard 
action. Interestingly, the topological lattice action approaches the continuum limit of the 
step scaling function from below, while the standard action approaches it from above. 

In this paper, we combine the standard and the topological action to a non-topological 
constrained action. The relative angle between neighboring spins is again limited by a 
maximal angle 5, but allowed configurations are now assigned the standard action value. 
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Similar actions with a constraint have been used before in various contexts [21-34]. Here we 
optimize the constraint angle 5 to reduce the cutoff effects in 2-d O(N) models. Remarkably, 
the cutoff effects of a variety of physical quantities including the step scaling function, the 
renormalized coupling g-^, and the mass in the isotensor channel, as well as the vacuum 
angle ^-dependence of the mass gap for N = 3, turn out to be at most a few per mille, 
even on rather coarse lattices. This provides us with a very simple nearest-neighbor action 
that can be simulated very efficiently with the Wolff cluster algorithm [10, 11]. In fact, 
the optimized constrained action reduces cutoff effects at least to the same extent as the 
parametrized classically perfect action, but is a lot simpler. A mixed action, which combines 
the standard quadratic with a large quartic term suppresses cutoff effects equally well. 

The paper is organized as follows. In section 2 we introduce the various actions to be 
considered in this work. Section 3 contains the investigation of the corresponding cutoff 
effects in the 2-d 0(3) model first at vacuum angle 9 = 0, and then also at 9 ^ 0. In 
section 4, we study the 2-d 0(4) and 0(8) models in a similar manner, concentrating 
on the step-scaling functions. Section 5 addresses the cutoff effects in the N = oo limit 
analytically. Finally, Section 6 contains our conclusions. An analytic investigation of cutoff 
effects in the 1-d 0(3) model is relegated to appendix A, while some technical details of 
the N = oo calculation are presented in appendix B. 

2 Lattice Actions for O(JV) Models 

In this section we introduce various actions for O(N) models. In later sections we will 
compare the corresponding cutoff effects in order to identify a highly optimized lattice 
action. In the continuum, the action of the 2-d O(N) model is given by 



where e(x) = (e±(x), e2(x), ■ ■ ■ ,ejy(x)) is an A r -component unit-vector field, and g is the 
dimensionless coupling constant. Just as non-Abelian gauge theories in 4 space-time dimen- 
sions, 2-d 0(N) models (with N > 2) are asymptotically free and have a non-perturbatively 
generated mass gap. For N = 3 one can define the topological charge 



which is an element of the homotopy group n^S* 2 ] = Z. In that case, one can introduce a 
vacuum angle 9 and add i9Q[e\ to the Euclidean action. Interestingly, as we have recently 
demonstrated, 9 is a relevant parameter and there are distinct continuum theories for each 
value 9 G [0, n] [35]. This conclusion has been further supported by [36]. For N > 3, on 
the other hand, 2-d O(N) models are topologically trivial. 

The standard action for the 0(N) model on a 2-d square space-time lattice takes the 
form 




(2.1) 




(2.2) 




(2.3) 
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where e x is an iV-component unit-vector associated with the lattice site x, and ft points 
from x to the neighboring site x + fi in the ^-direction. In the classical continuum limit 
the standard action with (3 = 1/g 2 reduces to the continuum action of eq. (2.1). According 
to Symanzik's effective theory, the standard action is expected to have cutoff effects of 
0{a 2 ). Recently, it has been observed that large logarithmic corrections, which also result 
from Symanzik's theory [18], can mimic the apparent 0{a) behavior that was observed in 
numerical simulations [16, 17]. 

Recently, we have performed detailed investigations of so-called topological lattice 
actions [20], which are invariant against small deformations of the lattice field. Here we 
investigate a topological action that constrains the angle between neighboring spins by 
a maximal value 5. The topological action vanishes, i.e. 5 top [e] = 0, if the constraint is 
satisfied for all nearest-neighbor spin pairs, i.e. e x ■ e x+ f l > cos 5, and is infinite otherwise. 
This action has already been considered earlier in [22-24, 26, 30], without emphasizing its 
topological features. Since the topological action vanishes for all allowed configurations, it 
has no meaningful classical continuum limit and cannot be treated in perturbation theory. 
Still, when one sends <5 — ?• 0, one approaches the correct quantum continuum limit [20, 
26], which demonstrates that universality does not rely on the classical continuum limit. 
Interestingly, for the topological action the sign of the cutoff effects of some observables is 
opposite to the one resulting from the standard action. 

In [35] we have used this observation to construct an optimized action with extremely 
small cutoff effects. The resulting constrained action combines the standard and the topo- 
logical action such that 

S con [e\ = ^2s(e x ,e x+fi ), (2.4) 

X,fl 

with 

s(e x ,e x+fl ) = /3(1 - e x ■ e x+ p) for e x ■ e x+ f t > cos 5, (2.5) 

and s(e x ,e x +ix) = oo otherwise. For 5 = ir the constrained action reduces to the standard 
action, while for (3 = it turns into the topological action. By optimizing the constraint 
angle to the value cos 5 = —0.345, we have been able to reduce the cutoff effects below the 
per mille level in an investigation of lattice #-vacua in the 2-d 0(3) model [35]. When 5 takes 
the optimized value, the constrained action -Sconle] turns into the optimized constrained 
action S oca [e\. In this paper, we construct such actions also for N ^ 3. 
Finally, we consider the following quadratic plus quartic mixed action 

Smixte] = ^2 ~ ^ x ' + T( 1 ~ ^ ' ^+a) 2 ] • ( 2 - 6 ) 

X,fl 

For given /3, we will adjust 7 to reduce the cutoff effects, thus optimizing the action. When 
7 > (3 2 the mixed action can again not be treated in perturbation theory. When 7 takes the 
optimized value, the mixed action 5 m i x [e] turns into the optimized mixed action 5 oma [e|. 

Berg and Liischer have introduced a geometric definition of the lattice topological 
charge [39]. In this definition, each lattice plaquette is divided into two triangles, as 
illustrated in figure 1. The spins e x , e y , and e z at the three corners of a lattice triangle t xyz 
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Figure 1. Decomposition of the square lattice into triangles t xyz . The topological charge density 
A xyz /4n of the 2-d lattice 0(3) model is given by the oriented area A xyz of the spherical triangle 
defined by the three spins e x , e y , and e z at the corners of t xyz . 

define the corners of a spherical triangle on S 2 . The oriented area A xyz of the spherical 
triangle is given by 

A xyz = 2ip G (— 2ir, 2tt], X + iY = r exp(z^), 

X = 1 + e x ■ e y + e y ■ e z + e z ■ e x , Y = e x ■ (e y x e z ). (2.7) 

The lattice topological charge is defined as the sum of the oriented areas A xyz over all 
triangles t xyz , normalized by the total area 4-7T of the sphere S 2 , i.e. 

txyz 

The decomposition of the square lattice into triangles illustrated in figure 1 is invariant 
under 90 degrees rotations and under translations by an even number of lattice spacings. 
When a configuration is shifted by just one lattice spacing, the lattice plaquettes are divided 
into two triangles in a different manner than before, and thus the topological charge may 
change. However, this does not happen for sufficiently smooth configurations. In particular, 
one can show that a nearest-neighbor constraint angle 5 < ir/2 (as used in the action S COQ ) 
leads to a completely translation invariant topological charge. In that case, just as in the 
continuum, different topological sectors are separated by infinite-action barriers. 

3 Numerical Study of Cutoff Effects in the 2-d 0(3) Model 

In this section we investigate the cutoff effects of a variety of physical quantities. In par- 
ticular, we address the question to what extent an action that was optimized to reproduce 
the continuum limit of a particular quantity automatically improves the scaling behavior 
of other observables. 
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3.1 The step scaling function at 9 = 

The dimensionless physical quantity 

u = Lm(L) = L/£(L), (3.1) 

is defined as the ratio of the spatial size L and the finite- volume correlation length £(L) = 
l/m(L). Based on this, one defines the step scaling function (with scale factor s) [9] 

a(s,u) = sLm(sL), (3.2) 

which is also known analytically [5]. By measuring the mass gaps m(L) and m(sL) in a 
Monte Carlo simulation on a lattice with L/a points in the periodic spatial direction, one 
obtains the lattice step scaling function T,(s,u,a/L) = sLm(sL,a/L). 

Figure 2 compares the cutoff effects of T,(2,uo,a/L) at uq = 1.0595 for five different 
lattice actions. While the standard action and the actions .0(1/3) and D(— 1/4) [40] (which 
contain additional diagonal neighbor couplings ), approach the continuum limit from above, 
the topological action approaches it from below. By optimizing the constraint angle 5 such 
that E(2,iio,a/L = 1/10) assumes its continuum value a(2, uo) = 1.26121035, one obtains 
cos 5 = —0.345. The resulting optimized constrained action has extremely small cutoff 
effects in the per mille range also at other (not too coarse) values of the lattice spacing. 
Based on the analytic results of [40] the lattice data have been fitted to the expression 

a 2 

Z(2, u, a/ L) = a(2,u) + [B log 3 (L/a) + Clog 2 (L/a) + .. .) , (3.3) 

which gives good agreement with the analytically known continuum result. The above 
analytic form is justified theoretically for the standard action as well as for the actions 
.D(l/3) and D(— 1/4). Although, strictly speaking, it may no longer be justified to use 
the expression of eq. (3.3) to fit the topological action data, one again obtains excellent 
agreement with the continuum result. Interestingly, in the range of correlation lengths 
considered here, the lattice artifacts of the topological action are smaller than those of 
the standard action. In fact, for the standard action at L/a = 64 the sub-leading term 
proportional to log 2 (L/a) is still larger than the leading term proportional to log 3 (L/a), 
while this is not the case for the topological action [20] . Based on the fitted values of B and 
C, the lattice artifacts of the standard action are smaller than the ones of the topological 
action only for correlation lengths larger than about 5 x 10 4 a. Without re-adjusting 5, using 
the optimized constrained action (which was optimized for no = 1.0595 at L/a = 10), the 
cutoff effects of the step scaling function are extremely small also for other values of u. 
Of course, one could also re-adjust S for each value of L/a (always optimizing at uq). 
Interestingly, the optimal value of 6 is rather insensitive to L /a (as long as L/a is not too 
small) . 

Let us also pursue this more elaborate alternative optimization strategy, however, now 
applied to the mixed action. For each value of /3, the parameter 7 of the mixed action has 
been optimized to reproduce the analytically known continuum limit of the step 2 scaling 
function, i.e. T,(2,uq, a/L) = a(2,uo) = 1.26121035. The optimal values of (3 and 7 are 
shown in figure 3 for various lattice sizes L/a. They are also listed in table 1. 



-6- 



0.04 



3 

w 

to 




Figure 2. Cutoff dependence of the step scaling function E(2, ito, a/L) for five different lattice 
actions: the standard action, the actions D(l/3) and D{— 1/4) wii/i additional diagonal neighbor 
couplings, the topological lattice action of [20], and the optimized constrained action with cos 5 = 
—0.345. The lines for the standard action are fits based on eq. (3.3) . The horizontal line represents 
the analytic continuum result of [5]. 
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10 


0.31664 


1.18790 


1.059500(6) 


1.261218(12) 


0.000008(16) 


16 


0.51378 


1.13480 


1.059506(9) 


1.261204(12) 


-0.000016(19) 


32 


0.73962 


1.10276 


1.059505(9) 


1.261167(17) 


-0.000011(22) 


64 


0.95450 


1.06820 


1.059501(15) 


1.261205(36) 


-0.000007(44) 



Table 1. Optimized parameters f3 and 7 for the mixed action such that 2(2, uq, a/L) = cr(2, uq) = 
1.26121035. The resulting value of u is very close to the desired uq = 1.0595, and E(2,ii, a/ L) devi- 
ates very little from the continuum value <r(2, uq). The last column is the deviation of S(2, uq, a/L) 
(obtained by extrapolation from u to Uq) from the analytic result. 



Using the couplings of table 1, figure 4 shows the resulting cutoff effects of the step 
3 scaling function T,(3,UQ,a/L) compared to the analytically known continuum result 
<r(3, no) = 1.439574. Without any further adjustable parameters, the mixed action op- 
timized in this way automatically leads to a drastic reduction of the cutoff effects of 
T,(3,u ,a/L). 

3.2 The step scaling function at 6 7^ 

While analytic results from the exact S-matrix exist only for 9 = and 9 = ir, the step 
scaling function is well-defined for all values of the vacuum angle 9. In [35] the step 2 
scaling function has been determined at 9 = 0,tt/2, and tt. In this way, the conjectured 
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Figure 3. Optimized parameters (3 and 7 for the mixed action such that S(2, uq, a/L) = <t(2, uq) = 
1.26121035. The dashed lines represent the tangents to lines of constant E(2, uq, a/L), which are 
obtained by measuring the derivatives du/d/3 and du/dj. 
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Figure 4. Cutoff dependence of the step 3 scaling function Y<(3,uo,a/L) for the standard action 
and the optimized mixed action. 



exact S-matrix has been indirectly verified with per mille level accuracy. It also has been 
shown that the step scaling function at 9 = ir/2 differs from the one at 6 = and tt. 
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Figure 5. Cutoff dependence of the step scaling function £(#, 2, uq, a/L) (with Lm(9,L) = uq = 
1.0595 ^ for the standard and for the optimized constrained action with cos<5 = —0.345, at 9 = tt/2 
and 9 = tt. The lines for the standard action are fits based on eq. (3.3). The horizontal lines 
represent the analytic result of [41] at 9 = tt, and the fitted continuum value for 9 = tt/2. 

This shows that 9 does not get renormalized non-perturbatively and it indicates that all 
values of 9 £ [0, n] are associated with distinct theories in the continuum limit. In view 
of the potentially devastating cutoff effects caused by dislocations, this is a non-trivial 
and somewhat unexpected result. Note, however, that a non-trivial ^-dependent spectrum 
is predicted by form-factor perturbation theory [37]. The spectrum was investigated by 
lattice Monte Carlo methods [38] using an imaginary 9 and analytic continuation. 

Figure 5, which is taken from [35], shows that without any further adjustments of <5, the 
cutoff effects of the optimized constrained action are automatically drastically reduced also 
at 9 = tt/2 and tt. The analytic result is known [41] for 9 = tt to be ct(tt, 2, uq) = 1.231064 
for u Q = 1.0595. 

We have also considered Lm(9,L) at fixed Lm(0,L) = uq = 1.0595, which is yet 
another physical quantity. Figure 6 illustrates the cutoff effects of the mass gap at 9 = tt 
for the optimized constrained action (with cos5 = —0.345), the optimized mixed action 
and the standard action. Although the cutoff effects are only in the per mille range for 
the optimized actions, the approach to the continuum limit is non-uniform since the lattice 
results undershoot the exact continuum value before they ultimately approach it from 
below. 

We have also measured Lm(7r/2, L) at fixed Lm(0, L) = no = 1.0595 using the opti- 
mized mixed action with the couplings listed in table 1. As shown in figure 7, the cutoff 
effects are again drastically smaller than for the standard action. 
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Figure 6. Cutoff dependence of the mass gap ratio m{9 = iv,L)/m(0,L) at fixed Lm(0,L) = 
Uo = 1.0595 for the standard, constrained (with cosS = — 0.345J, and optimized mixed action. The 
horizontal line is the analytic result. 
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Figure 7. Cutoff dependence of Lm(0,L) at fixed lm(0, L) = uq = 1.0595, for the standard and for 
the optimized mixed action at 8 = tt/2. The horizontal line represents the fitted continuum value. 
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Figure 8. Cutoff effects of the mass of the isotensor state for the standard and optimized mixed 
action (with the parameters listed in table 1). 



3.3 The mass of the isotensor state 

Besides the state in the isovector channel, the exact S-matrix also provides analytic results 
for the finite- volume mass rri2{L) in the isotensor channel (7 = 2). Its approach to the 
continuum limit is shown in figure 8 for both the optimized mixed and the standard ac- 
tion. Again, the optimized mixed action reaches the continuum limit much faster than the 
standard action. 

3.4 The renormalized coupling gn 

The renormalized coupling defined in terms of a truncated 4-point function at zero momen- 
tum (see eq. (5.44)), has been calculated from the exact S-matrix in [42] with the result 
<?r = 6.770(17). As illustrated in figure 9, without any further tuning of the optimized 
actions, the continuum limit of gn is approached much more quickly with the optimized 
constrained and the optimized mixed action than with the standard or the topological 
action. The results from the optimized actions suggest that the theoretical error given in 
[42] has been overestimated. Based on the last two points with the optimized mixed action 
here we obtain the value $r = 6.769(2). 

4 Numerical Study of Cutoff Effects in the 2-d 0(4) and 0(8) Models 

In order to investigate whether highly optimized local actions can also be constructed 
successfully for larger values of N, in this section we study the 2-d 0(4) and the 2-d 0(8) 
model. In these cases we limit ourselves to the optimized constrained action compared to 
the standard action. 
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Figure 9. Cutoff effects of the renormalized coupling for four different lattice actions. The 
thick horizontal line is the result = 6.770(17) obtained from the exact S-matrix (whose estimated 
theoretical error corresponds to the dotted horizontal lines) [4-2]. 

4.1 The step scaling function of the 2-d 0(4) model 

Figure 10 compares the cutoff effects of the step 2 scaling function for the standard action 
with those of the optimized constrained action. The constraint angle 5 has been optimized 
by demanding E(2, no, a/L = 1/10) = er(2, no), in this case for uq = 1, which yields cos 5 = 
—0.096. Again, the cutoff effects are drastically reduced in comparison to the standard 
action. Still, there are remaining tiny cutoff effects in the per mille range. The cutoff effect 
is not monotonic, after a small increase it starts to diminish only on an L/a = 32 lattice. 
The data for the standard action are taken from [18]. 

4.2 The step scaling function of the 2-d 0(8) model 

Figure 11 shows the analogous results for the 0(8) case. In this case, the constraint is 
cos 5 = 0.217, again for no = 1.0595. The data for the standard action are from [17]. 

5 Analytic Study of Cutoff Effects in the 2-d O(N) Model at N = oo 

The step scaling function for the 0(N) non- linear sigma model in the limit N — > oo has 
been studied in [17, 43, 44]. Here we generalize these calculations for the case of the 
constrained action and of the mixed action. Note that here we use lattice units, "a = 1" . 
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Figure 10. Cutoff effects of the step-scaling function in the 2-d 0(4) model. The horizontal line 
indicates the exact result in the continuum limit, er(2, 1.0) = 1.208658. The optimized constrained 
action has cutoff effects in the fraction of a per mille range. 
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Figure 11. Cutoff effects of the step-scaling function in the 2-d 0(8) model. The horizontal line 
indicates the exact result in the continuum limit, a (2, 1.0595) = 1.345757. The optimized constrained 
action has cutoff effects in the fraction of a per mille range. 
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5.1 Constrained action in the large iV limit 

The action and the partition function in the constrained case are given by 

5 [ e 1 = ^E(^) 2 , (5-1) 

Z= /exp{-5[el}n^1-l)n ( e -(^) 2 )- (5-2) 

e X X,fJ, 

Here is the step function, denotes the forward lattice derivative, and J g denotes 
/ Tlx d<4- As usual, the N — > oo limit is taken with (3/N = l/f fixed. 

For the standard action the distribution of (d^e x ) 2 in the large iV limit approaches a 
(5-function at (c^e^) 2 = f jd (where d is the number of Euclidean dimensions), as can be 
easily obtained from perturbation theory. On the other hand, for the topological action 
(with only the G-constraint, i.e. with / = oo) it goes also to a 5-function at (dfj,e x ) 2 = e. 
One expects that for e > f/d the constraint is irrelevant (it is above the narrow peak) 
and / alone determines the physics in the large N limit. On the contrary, at e < f/d the 
constraint dominates the physics. It is also clear from these considerations that at N = oo 
the 0(e — (d^e x ) 2 ) factors can be replaced by a strict constraint 5(e — (c^e^) 2 ). 

Introducing the auxiliary variables a x , rj x ^ one gets 

Z= ( exp{-5[el+i(a,e e -l)+i(r ? ,e-(9e) 2 ))^-. (5.3) 

Je,a,ri 1 J 7? — zO 

Here the auxiliary variables are integrated over the real axis. 

Note that e = 2(1 — cos 5) where 5 is the maximal allowed angle between the near- 
est neighbors. The allowed range is < e < 4. For e > 4 one recovers the standard 
(unconstrained) action. 

Rescaling a and r\ as 

a x -» 2 Nax ' ^ ~^ \ Nr lxti , (5.4) 

and then shifting the integration contour 1 

a x ih + -yL , rj Xfl -> -iVf, + -^L , (5.5) 
v TV viV 



one obtains a form suitable for the 1/N expansion: 
N 



S eS = - 



E (7 + v ») (^) 2 + h £ e 2 x - Vh - eV ]T v 
-iylE^-^+E^-Wlj+j;^^*!) , (5-6) 



1 T)\xe to eq. (5.3) the integration line of rjxy, can be shifted only downwards, i.e. one should have > 0. 
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Note that in eq. (5.5) we have allowed for different values of in different directions. 
In symmetric volumes V = L d they can be replaced by a single v. For the step scaling 
function one should take the strip geometry, oo x L^ 1 . Here one should expect different 
«u's in the time- and spatial directions. 2 

3 



After integrating out the spin fields we get" 



tr 



\nS -Vh-Ve^v^ 



+ O (Vn) , (5.7) 



with 

Sxy = h5 X y + ^ \25 X y — 5y^ X +jX — Sy^ — p] , (^-8) 



f 

7 

For the inverse we have 



Wp = - + . (5.9) 



D(p) = Y / w^(p) + h, (5.11) 

where the sum is over p^ = 27rn M /L /J , = 0, 1, . . . , — 1 and where K^(p) = 2(1 — 
cospp). 

Note that the l/(r] — i0) denominator (i.e. the logarithmic term in eq. (5.6)) does not 
enter in the N = oo limit, i.e. the O-constraint is effectively replaced by a ^-constraint, as 
discussed above. 

Eq. (5.7) in leading order yields the gap equations 



1 1 , 



V ^ D(p) 

Consistency of these equations gives the relation 

h + eJ2 w v = 1 - ( 5 - 14 ) 

The finite volume mass gap m(L) in the strip geometry oo x L d ~ l in leading order 1/N 
is given by 

m = 2sinh^^, (5.15) 



2 Note that one could avoid this complication by modifying the constraint to an "isotropic" one, 
0^e— (de x ) 2 ^j, meaning that the average of (d^Sx) 2 over the 2d nearest neighbors should not exceed 
e. 

3 In the leading order in 1/N we will not need the terms depending on the auxiliary fields a x , rj xfl . Hence 
they are not written out explicitly here. 
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where 

Defining 
one has 



h 



m z = — . (5.16) 



u = m(L)L, (5.17) 



u u I u 2 



m = m (u,L) = 2 sinh — = _M + _ + ...] . (5.18) 



The continuum limit is approached by L — > oo keeping u fixed. 
5.1.1 The isotropic case 

In a symmetric volume L d (or with the modified "isotropic action") one has i>« = i> and 
the gap equations (5.12), (5.13) yield 

V ^ K{p) + m§ ~ m§ + ed ' (5-19) 



as well as the relation 



v + j)ed + h = l. (5.20) 

Since one should have h and v non-negative, this relation has a solution only for e < f /d. 
As discussed above, for larger e the kinetic term is the relevant one and the constraint can 
be omitted. As expected, the gap equation (5.19) has a solution only for e < 4 since for 
the standard action < K(p) < 4. 

For the standard (unconstrained) action the gap equation reads 



-T. 



V^K(p) + mi f 



(5.21) 



In the infinite volume eq. (5.19) gives 

d d p 



(271-)^ K(p) +ml ml + ed 



(5.22) 



From now on we consider the d = 2 case. In this case the integral on the l.h.s. 
of eq. (5.22) diverges logarithmically like l/(4-7r) log(a 2 mo), where we restored the lattice 
spacing. For the standard action this gives the well known result 

const n-it , . 

m w e~ 2lT/f . (5.23 

a 

Since ttiq oc exp(— 2ir/e) in the continuum limit e — > the rriQ term can be neglected 
on the r.h.s. of eqs. (5.19), (5.22), and in the continuum limit one recovers the universal 
result. The cutoff effects will be, however, slightly different. To determine them for the 
step scaling function one can solve the gap equation (5.19) numerically in a finite volume 
for different lattice spacings (i.e. different L/a). It turns out that the leading 0(a 2 ) cutoff 
effects can be calculated analytically. 
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5.1.2 The lattice artifacts at AT" = oo for the step scaling function 

Consider the d = 2 case in the strip geometry oo x L and introduce 

P=^, (5-24) 

00 = P K x {p) + ml . (5.25) 

Writing eqs. (5.12), (5.13) for Lt = oo after integrating out po one obtains a set of two 
equations 

L p v^+4) 

I^^lM= = e „ , (5.27) 

where from eq. (5.14) one has wq = l/[mg + (p + l)e]. 

We are interested in the approach of the step scaling function u'(u, L) = m(2L)2L for 
fixed u(L) = m(L)L to the continuum limit a/L — > 0. Analytic results are obtained using 
the results of Caracciolo and Pelissetto [44]; technical details are deferred to appendix B.l. 
The result is of the form 

u' = u' OQ + j^v{u,z) + ... , (5.28) 

where the continuum limit u'^ = <r(2, u) is (the same for all actions) given by the solution 
of the relation 

/o(n) = / (n' oo ) + ^ln2, (5.29) 

where the function fo(u) is specified in eq. (B.12). The lattice artifacts depend on the 
lattice actions; the leading artifacts are specified by the function u(u, z) where, as will be 
seen later, it is convenient to express the InL dependence through z defined by 

z = z(u,L) =/ (u) + ^-lnL. (5.30) 

Z7T 

Considering first the standard action, the leading lattice artifacts are linear in z: 

Vstd(u,z) = t (u) +t±(u)z, (5.31) 

where the functions to(u),ti(u) are specified in eqs. (B.22), (B.23). The coefficient t\(u) 
is positive, so the asymptotic approach to the continuum limit is always from above. Nu- 
merical evaluation shows that to(u) is negative for all u, hence the two terms compete with 
each other. It turns out that v s td(u,L) is positive at u > 0.6357 for any L. Below this 
value the approach for reasonable L seems (misleadingly) to be from below. For small u 
the asymptotic behavior sets in only at very large L/a, e.g. at u = 0.3 the coefficient 
v s t&{u, L) becomes positive only for L/a > 10 5 . 

Turning to the constrained action, one can calculate the infinite volume correlation 
length £ for this case. For large L one gets from eqs. (B.3), (B.5) 

1 = Z + 0{1/L 2 ) . (5.32) 
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1/N 

Figure 12. u(e, L) for the O(N) topological action at large N values. The line is a cubic fit to the 
5 data points. The two points at 1/N = show the results obtained from solving the gap equations 
numerically. The lower point is obtained from eqs. (5.26), (5.27). For comparison we also give the 
result for the "isotropic action", eq. (5.19), upper point. 

Using the asymptotic expression for large u = m(L)L ~ ML 

/o(u) = -^lnu + £ln2 + 0(e-") , (5.33) 

one obtains 

l = ^ln(x/32£)[l + <?(r 2 )] • (5.34) 

This is the same relation as for the standard action [43] when the coupling / is replaced 
by 2e. 

First, one can solve the coupled equations determining the step scaling function for the 
constrained action numerically as indicated in appendix B.l.l. We checked the numerical 
solution by direct Monte Carlo measurements for the constrained topological action (/ = 
oo) at e = 0.2, L = 10 for N = 10,20,30,40,60. The results are shown in figure 12. The 
value extrapolated to TV = oo is 0.225(1). It agrees with the solution of the coupled gap 
equations, u(0.2, 10) = 0.225410. For the anisotropy parameter in this case one obtains 
p = 0.950485. (Note that for the isotropic model, eq. (5.19) one gets a slightly different 
value, u(0.2, 10) = 0.230919.) 

One can also proceed analytically (see appendix B.1.2) and one finds that the lattice 
artifact u con (u,L) is now quadratic in z: 

v CO n(u, z) = io(u) + h(u)z + i 2 (u)z 2 , (5.35) 
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where the functions ti{u) are given in eq. (B.29). Comparing this to the result for the 
standard action, one has 

f 2 (u) = -8ti(u), (5.36) 

i.e. the asymptotic approach of u' to its asymptotic value for the constrained action is 
always opposite to that of the standard action, so that the approach in this case is always 
from below. This behavior is essentially the same as observed for the 0(3) case. 

We remark that the In 2 L behavior found above is not in contradiction to the analysis 
of ref. [40], since the constrained action does not belong to the class of actions considered 
there. 

5.2 Mixed action in the large N limit 

Here we consider the mixed action 

Sri*® = \ EO 9 /^) 2 + \ EPm4) 2 ] 2 • (5.37) 

X,fl X,fl 

It has the advantage over the constrained action that the cutoff effects can be changed 
continuously, and in particular, tuned to zero for a given quantity, similar to the 0(3) case 
discussed in section 3. 

We shall take the large N limit as 

N 2N 
P = -> 7 = — • (5-38) 



The 1/N expansion is obtained similarly to the case of the constrained action and is 
described in appendix B.2. 

Introducing the effective coupling / = /(/, k) by 



11 / 1 1 

7 = 27 + Vip + ^' (5 ' 39) 



similarly to the constrained case (cf. eqs. (5.32)-(5.34)) one finds for the infinite volume 
correlation length 

i = ^ln(V32 0[l + O(r 2 )] • (5.40) 



Besides / we introduce 



and q = - (r + \J r 2 + 4^ . (5.41) 



H 
1 



The fixed points in the coupling space are those where /(/, k) = 0, i.e. the boundaries 
of the first quadrant of the /, K plane. The continuum limit of the step scaling function, 
ti'oo = a(2,u) is of course universal, but the cutoff effects depend on the ratio r, in general 
on the particular path r(/) along which one approaches to the continuum limit. 

In appendix B.2 it is shown that the leading artifact is of the form 

^mix(«, z) = T (u) + T x {u)z + T 2 (u)z 2 , (5.42) 
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Figure 13. Deviation of 11(2, u, a/L) from the exact value for the standard and quartic action. 



with functions Ti(u) specified in eqs. (B.53), (B.54). 

One notes that the coefficients appearing in Tq(u) and T2(u) depend on the ratio r. 
For r = constant the cutoff coefficient v m ix(u, z) is in general a second order polynomial in 
z (i.e. in InL, cf. eq. (5.30)). The exception is the standard action (r = oo) where it is a 
first order polynomial. In this case one recovers eqs. (5.31) and (B.23), as expected. The 
purely quartic action is obtained by setting r = 0. 

The leading cutoff effect is positive for the standard action, and (for fixed r) negative 
in other cases. Figure 13 shows the deviation 11(2, u, a/L) — a(2,u) for the standard and 
quartic action (with k = oo and / = oo, respectively.) 

As one can see from eqs. (5.41), (B.54), for large z one can cancel the artifact by 
choosing r a function of / so that asymptotically q 2 = 8z + 0(1). This gives for large L 
the optimal path 

(5.43) 

Note that the cancellation between the z and z 2 terms takes place for all values of 
u and scale factors s, hence the lattice artifact T,(s,u,a/L) — a(s,u) is expected to be 
significantly reduced for all s and u once it is removed, say, for s = 2, u = 1. 

The optimal path in the /, k plane shown in figure 14 is determined from the condition 
T,(2,u,a/L) = a(2,u) using eqs. (B.43)-(B.45). The dashed line corresponds to n 2 /(2f) = 
4. Figure 15 shows the deviation 11(3, u, a/L) — a(3,u) for the optimized values /, k. 
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Figure 14. The curve of the optimal couplings determined from E(2, u, a/L) at u = 1 for JV = oo. 
The solid curve shows the results for L/a = 10, ... , 2000, the dashed line shows the asymptotic 
dependence, K = \/8f. 
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Figure 15. Deviation of E(3, u, a/L) from the exact value for the standard, quartic and optimized 
mixed action. For the latter the couplings are determined to have E(2,ii, a/L) = a(2,u). 
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5.2.1 The renormalized 4-point coupling 

Here we consider the renormalized 4-point coupling at zero momentum, gn, defined through 
the Binder cumulant: 

VM 2 ^ ^ _ ^ 

x,y,x',y' 

£ = • 4) , (5.45) 

and M is some renormalized mass: 

M = m + 0(1/N) . (5.46) 

Since g-& = 0(1/N) for large N one needs a calculation at order 1/iV which is described in 
appendix B.2.1. 

Taking the infinite volume limit and the continuum limit, the lattice artifacts are given 
by (see appendix B.2.2) 4 

Ng R (aM) = 8vr + a 2 M 2 A{q, C) + C(a 4 M 4 ) , (5.47) 

where C = - ln(M 2 /32) and 

A(q,C) = 27T + -^- (k + —^)c + -^C 2 . (5.48) 
1 + q z \ 1 + q z J 1 + q z 

Here q is given by eq. (5.41). 

We see that the artifacts in gn behave very similarly to the step scaling function 
artifacts (except that they are of opposite sign) in that the leading artifact is 0(a 2 M 2 ) 
times a second order polynomial in ln(aM). 

Approaching the continuum limit along the curve K?/(2f) = 4 (cf. eq. (5.43)), for 
which the logarithmic terms in the step scaling function cancel, we have q 2 = 1 + 8w, and 
a similar cancellation occurs here. 

The numerical evaluation of Ngu in the infinite volume limit is described in ap- 
pendix B.2.3. Figure 16 shows the deviation Ng^(aM) — 8ir as a function of aM. 



6 Conclusions 

We have investigated a variety of lattice actions in numerical simulations of 2-d 0(3), 0(4), 
and 0(8) models, as well as analytically at iV = oo. By optimizing the constraint angle 5 
or the quartic coupling 7 in very simple nearest-neighbor actions that suppress large field 
fluctuations, we have obtained a class of lattice actions with extremely small cutoff effects, 
often in the fraction of a per mille range. The simplicity of these actions makes them very 
useful for numerical simulations. In particular, the Wolff cluster algorithm is applicable in 
a straightforward manner. 

4 Here we re-introduce the lattice spacing a. 
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It is pointed out in [20] and also observed here, that in order to belong to the corre- 
sponding universality class an action is not required to have the standard classical contin- 
uum limit. 

Interestingly, both the constrained action and the mixed action have unexpected cutoff 
effects at N = oo. In particular, the leading artifact in both cases is ~ a 2 log 2 (aM) while 
for the class of lattice actions considered in [40] in the framework of Symanzik's effective 
theory the leading artifact is ~ a 2 log(aM). 

Although this might appear as a puzzle at first sight, in fact Symanzik's effective 
theory approach is valid also for the actions considered here. In this framework the lattice 
artifacts are described by the effective continuum Lagrangian 

£ arti = a 2 C0 4 (6.1) 

where the c- number coefficient C is coupling-dependent and O4 is a local 0(iV)-invariant 
operator. (In fact the complete artifact Lagrangian is a linear combination of several terms 
of this form.) Approaching the continuum limit along a given curve in coupling space, C 
becomes a function of the inverse correlation length aM. On the other hand the matrix 
elements of the operator O4 are cutoff-independent but depend on the physical parameters 
(u for the case of the step-scaling function). We know how to calculate the asymptotic 
form of the dependence on aM only if coupling constant perturbation theory is applicable. 

The actions considered in [40] were perturbative — close to the continuum limit the 
fluctuation of the corresponding field was determined mainly by the quadratic part of 
the action. On the other hand, for the constrained action (for sufficiently small e) these 
fluctuations are restricted by the constraint, while for the mixed action (for sufficiently 
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large 7) by the quartic term coming from 7(c^e) 4 , hence the situation is non-perturbative. 
(This is obvious for the topological action, which is zero for the allowed configurations.) 

The large N results suggest that the coefficient C receives non-perturbative contribu- 
tions, which behave as oc log 2 (aM) for N ->■ 00. This scenario is supported by the fact that 
for these actions the log 2 (aM) and the log(aM) terms have the same li-dependence (namely 
u 2 , cf. eqs. (B.28), (B.54)) as the leading term for the perturbative actions (coming from 
the same operator). This issue deserves further investigation. 

At large N for the standard action the squared fluctuation of the relative angle is 
roughly given by /, while in the case of the mixed action for the purely quartic case 
(/3 = 0, i.e. / = 00) it is given by k. If by taking the continuum limit the ratio r = n/f 
approaches infinity, then the fluctuations are dominated by the quadratic part of the action. 
In this case one observes indeed an a 2 log(a) artifact, as for the other perturbative actions. 
In particular, this is the case when one moves along the optimal curve k = y/8f, where 
even the leading a 2 log(a) artifact cancels. By taking the continuum limit along a straight 
line, r = const both the quadratic and the quartic terms are relevant. 5 This case cannot 
be treated in perturbation theory, and produces a leading a 2 log 2 (a) artifact. 

It is remarkable that a very simple nearest-neighbor action can reduce cutoff effects 
that drastically, and it is obvious to ask whether this success extends to other interesting 
asymptotically free lattice field theories, including 2-d CP(N — 1) models, 4-d Yang-Mills 
theories, or even QCD. Optimized actions are straightforward to construct in all these cases, 
and investigations in this direction are currently in progress. Of course, there is no exact 
result that helps finding the optimalized values of the parameters in these cases. However, 
as was found here, determining the optimal choice of the parameters by calculating one 
physical quantity very precisely leads to an action where lattice artifacts are small for a 
large class of other physical quantities as well. 
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A Analytic Study of Cutoff Effects in the 1-d 0(3) Model 

For the constrained action in one dimension the eigenvalues of the transfer matrix (with a 
properly chosen normalization factor) are given by 

K(P,c)= [ dzJl'-Vp n (z) (A.l) 



The standard action corresponds to c = cos 5 = — 1. 
With the notation x = 1/(3, e = exp((c — l)/3): 



A = x - ex , (A.2) 



Ai = x — x 2 — e(cx — x 2 ) , (A. 3) 



A 2 = x - 3x 2 + 3x 3 - e ( ^(3c 2 - l)x - 3cx 2 + 3x 3 ) , 



A 3 = x - 6x 2 + 15x 3 - 15x 4 



( ^(5c 3 - 3c)x - ^(15c 2 - 3)x 2 + 15cx 3 - 15x 4 ] 



The excitation energies of the rotator (in lattice units) are 



aE n = -log (^\. (A.6) 

The lattice spacing is determined here by aE\ = 1. 
The ratios of the energies are given by 

r n = ^. (A.7) 

In the continuum limit (c fixed, (3 — > oo) they approach the well-known result for the 
quantum rotator, 

lim r n = -n(n + 1) . (A. 8) 

(3— >oo 2 



In particular 



r 2 = 3 - x 2 - -x 3 + . . . (A. 9) 

r 3 = 6-5x 2 -yx 3 + ... (A.10) 



In the case of the topological action, /3 = 0, c — > 1, one has 



r 2 = 3 + i(l-c) + ^l-c) 2 + ... (A.ll) 
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Figure 17. The ratio of the mass gaps for the 1 = 2 and 1=1 states as a function of the lattice 
spacing for the standard, topological, and constrained action. For the latter the constraint is fixed 
at E x a = 0.1 



5 25 

r 3 = 6 + -(l- C ) + -(l-c) 2 + ... (A.12) 

Comparing eqs. (A. 9), (A. 10), (A. 11), (A.12) one sees that the discretization errors for the 
standard action and for the the topological action have opposite signs. Fixing the constraint 
e.g. to c = 0.4 one obtains the resulting artifacts for T2 which lie in between those for the 
standard and the topological action (cf. figure 17) An interesting possibility is to choose 
c = c(/3) in such a way that the first ratio is exact, r2(/3,c(/3)) = 3. Of course, this 
procedure does not eliminate the discretization errors in other quantities, but it turns out 
that it improves the convergence to the continuum limit. Figure 18 shows r^{j3,c{(3)) — 6 
as a function of aE±, which in a log- log plot looks nearly linear. From the numerical values 
one finds that r^fi, c((3)) — 6 oc a a , with the power a ~ 2.77. 
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Figure 18. The deviation of the ratio of the mass gaps for the 1=3 and 1=1 states from its 
continuum value. The constraint c = c(/3) is chosen so that EijE\ = 3. 



B Technical Details of the N = oo Calculation 

B.l Analytic behavior of the step scaling function: constrained action 

Introducing 



L p V w ( w + 4 ) 

I 2 (p,u,L) = jJ2 , " (B-2) 
L p V w ( w + 4 ) 



(with mo given by eq. (5.18)) one obtains 



h(p,u,L)= , (B.3) 

m o + \P + 1 ) e 

,2 



mo + pe 
m l + {p + !) e 



J 2 (p^,L)= _ pL ^- . (B.4) 



The asymptotic expansion of I\, I 2 can be obtained from the results of Caracciolo and 
Pelissetto [44] (expressing the expansion in terms of u = mL instead of m^L used therein). 
One has asymptotic expansions 

h(p,u,L) =z+ jpB(u,z) + ... , (B.5) 

h(p,u,L) = 1 + -Ld(u,z) + ... , (B.6) 
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where we recall z is defined by 



z = z(u,L) = f (u) + — InL. 

Z7T 



and 



B(u, z) = fi(u) + -u 2 f {u) - ^u 2 z - y 
D(u, z) = f 2 (u) - -u 2 f (u) + ^u 2 z + |^ 
Here p\ appears in the expansion of p(u, L) 



z + ufo(u) + 



4"7T 



p(u,L) = l + j^ Pl (u,z) + 0(l/L 4 ) . 



(B.7) 

(B.8) 
(B.9) 

(B.10) 



Note that the isotropy is recovered (p — > 1) for L — > oo, as one can see from eq. (B.4) and 
from the infinite volume limit 



l2(p,u,oo) = — arctan (-Jp) ■ 

7T 

The remaining L dependence in p\ is assumed to be only logarithmic. 
The functions j%(u) are given by 

am = k + i 



1 

12 



f 2 (u) = 2nd 



IT U U . 

+ - + — (2k- 1), 



{-, 

\2ttJ 6 2 8tt 
k = je — In 7T H — In 2 . 



12 ~ lOTT 

2 



3 V27T 



We also need the derivative of fo(u): 



/<5(«) 



2u 2 (2vr) 3 L \2n 



Here 6 



G (x) = 

n=l 

oo 

Gi(x) = 



1 



y/n 2 + x 2 n 



n=l 

oo 



\J n 2 + 



„2n 



2n 



2\-3/2 



n=l 



(B.ll) 



(B.12) 
(B.13) 

(B.14) 
(B.15) 

(B.16) 

(B.17) 
(B.18) 
(B.19) 



3 Note there is an error in eq. (A.5) of [44]: it should read Gi(a) = 2J2T=i i FFr ( 2 fe te )C(2fc + 1) (|) 2fe+2 
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Note for further reference that /q(u) < 0. 

For the step scaling function we are interested in u'(u, L) = m(2L)2L for fixed u(L). 

Consider first the standard action. From the gap equation I\ (u,L) = 1/ f and I\ («', 2L) 
1/ f for u' = u'(u, L) one has the equation 



z + 



L 2 



fi(u) + \u 2 f {u) - \u 2 z 



z' + 



4L 2 



/l( </) + I^/ o( „')-^/V 



+ 0(1/L 4 ). (B.20) 



Expanding u' as in eq. (5.28) one sees that the continuum value u'^ is given by the solution 
to eq. (5.29). 

From here and eq. (5.30) one has for z' = z(u',2L) 



i' = z+^ti(u'Mu,z) + o(^J . 



(B.21) 



Because of this relation one can replace z' by z in the 1/L 2 terms in eq. (B.20). This is 
the reason for using z instead of InL. 

For the leading lattice artifact one obtains eq. (5.31) where 



t (u) 

h{u) 



l 



/0«o) 

i 



i 



i 



i 



fi(u) + -u 2 f (u) - - hiu 1 ^) + -u'iUu'^) 



2 ^ 12 

u -4^00 



(B.22) 
(B.23) 



Note that according to eq. (B.16) /q(u) is negative for all u. Since u'^ < 2u the 
coefficient t\{u) is positive. 



B.l.l Solving the coupled equations numerically 

Eliminating e from eqs. (B.3),(B.4) one gets (with uiq given by eq. (5.18)) 

(1 + p)I 2 (p,u,L) = mlh{p,u,L) + p. 
This determines p = p(u,L). Then one finds the corresponding e by 



1 



m, 



1 



h(p,u,L) 

Keeping e fixed one can calculate v! = u(e, 2L). 

B.1.2 Analytic form of the leading artifact 

Inserting the asymptotic expansions (B.5),(B.6),(B.10) we get 



(B.24) 



(B.25) 



2vr 
vr-2 



(2f 2 (u)-u 2 f (u)) , 



which is independent of z and 
1 Ph 



2e 2{I 2 - m 2 h 



($o(«) + $i(u)z + §2{u)z 2 ) + ... 



(B.26) 



(B.27) 
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where 



*o(«) = fi(u) + \u 2 f Q {u) - ipx(u) (uf Q (u) + 
1 



*i(«) = -g« , 

$ 2 («) = u 2 . 
One now obtains eq. (5.35) with 



(B.28) 



*(«) = ^ (**(«) ) ■ 



(B.29) 



B.2 1/iV expansion for the mixed action 

Here we discuss the mixed action given in eq. (5.37) in the large N limit given by eq. (5.38). 



1). 



The partition function is 

Z = ^exp | - Smbje]} Yl <5(4 
Introducing the auxiliary variables a x , f] x ^ one obtains 

Rescaling and shifting the integration contour as 



(B.30) 



(B.31) 



1 »t / .j . «x 
a x — )■ -iV ?/i H = 

2 I v 7 ^ 



(B.32) 



one gets 



^eff[e1 = -iV 



--2 

K v — \ 



X,fJ, 



(el - 1) 



. (B.33) 



Note that one has to allow for an anisotropy in because the lattice is not cubic, Lt ^ L s . 
After integrating out the spin fields we get an effective action in the auxiliary fields 



ScS = Y 



tr 



■2 

In R-Vh- V— 



4 A^/ "A* 
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with 



i 

Rxy = S X y + -—T X y i (B.35) 



S X y = kS X y + Wfj, [25 X y — Sy jX+ jj, ~ Sy^-fr] , (B.36) 



-Try — C^x^xy txy j (B.37) 
= [t/x/i ($sy — ^y,x+(i) + ^(x-fi)^ ($xy ~ <>y,x-fi)] > (B.38) 



/' 
1 



«V = j + v„ . (B.39) 

For the inverse S~y, and D(p) we have the previous expressions, eqs. (5. 10), (5. 11). 
In the leading order of 1/N expansion one obtains the following equations for h and 



v^wr 1 - (B40) 



This yields the relation 



-K 2 J2 V » W » + h = 1 - ( B - 42 ) 



2 

A* 



For the d = 2 case, with the same notations as in the constrained case (cf. appendix B.l) 
we get for L t = oo: 

h(p,u,L)=w , (B.43) 

I 2 (j>, u, L) = w (*jpvi + rn%\ , (B.44) 

and 

-K 2 (v + /9Vi)+mg = — . (B.45) 
2 Wo 

Eliminating vq, V\ from eqs. (B.43), (B.44) and (B.45) we can express the two couplings 
in terms of u, L and the remaining saddle point parameter, p: 

2 _ T 2 mA ^ 

~ ~ f, i „m- : — ZjTT 1 ! ' (BAb) 



k 2 (1 + p)I 2 - p - mg/i 
1 _ 2 J 2 - 1 

7 ~~ + ~2 ' — r 



~ = h + Z2 ■ ^7— ' ( B - 47 ) 



Inserting the asymptotic expansions (B.5), (B.6) into eq. (B.47) one obtains 

7 = 2 -i + G?)- (B48) 
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Solving this for z one finds 

z(u,L) = j + o(Jf) , (B.49) 

where the effective coupling /(/,«;) is defined in eq. (5.39). 

Inserting the asymptotic expansions 7 (B.5), (B.6), and (B.10) into eq. (B.46), and 
using eq. (B.49) one obtains 

2(2f 2 (u)-u 2 f (u)) 

Pl ^ = ! _ 2 , 2 ' ( B - 5 °) 

up to 0(L~ 2 ). Here q is a function of the coupling ratio r = k/ f given by eq. (5.41). Note 
that taking r = oo one should recover the result for the standard action, while taking r = 
that one for the pure quartic action. 

We now eliminate I 2 using eq. (B.46) and obtain from eq. (B.47) 

1 1 f. 9s 2 2m 2 ,) 

Inserting the asymptotic expansions one obtains 



1 1 1 

/ K 2 Z L 2 



fi(u) + \u 2 f (u) - \u 2 z 



;Pi{u) ( «/o(«) + ^ 



Using these relations one obtains the leading lattice artifacts given in eq. (5.42) with 
the functions Tj(u) given by 

Ti(u) = (^(u) - i*i(T4,)J , (B.53) 



$!(«) = -i« 2 , (B.54) 

o 

l + q z 

B.2.1 The renormalized 4-point coupling in leading order 

In order to compute at higher orders in the 1/iV expansion let us introduce a source in the 
action; 

S e g[e, J] = Ses[e\ + y~] J x ■ &x (B.55) 

X 

Then after integrating out the spin fields we get 

S cS [J] =S cS -^Yl UR~ l )*yJy ■ (B.56) 



7 Using eqs. (B.46) and (B.ll) one can show that p(u, 00) = 1. 
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The propagators of the auxiliary fields are given by the quadratic terms 

;,r (S- 1 TS- 1 T) + - 

4 



S eff , quad ratic = ~M (S^TS^T) + ^ , (B.57) 



where S X y is given by eq. (5.10). 
Defining the Fourier transforms 



v 

we have 8 



<*x = y'Z,j px 5(p), (B.58) 
v 

^ = ^E e * (<B+A/2) W. (B.59) 



1 Jp{x-y)+iqx r 

(TS ) xy = V iYl D(v) |- 5 (?) + (P + . ( R6 °) 



D(p) 

p,q 



so the quadratic term in the auxiliary fields is 



1 1 [q(g) - p^ri^q)] [a(-q) - p v f v ri v (-q)] , 

4V 2 f. D{p)D{r) ' P + q - [ 1 

Define {r = p + q): 

B ^-v^mwr (B62) 

W^Sjw. (B.63) 

Note all ii— functions are even in q. Then the leading order quadratic term can be written 

Seff.quad = ~^2[a(q)H{q)a(-q) - a(q)H ll (q)rj fl {-q) - a(-q)H fl (q)rj^ l (q) 
q 

1 K 2 

+ Vn(<i) H nv(<i)vv(-q)\ +^E^ 9 ^N- ( B - 65 ) 

Defining 

= S(g) - §^(?) . CB.66) 



Here p M = 2sin(p^/2). 
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we have diagonalized the quadratic part: 

1 1 



g 

+ Mi) { H ^i) ~ H " { H{ q ) {q) + l K%v ] U ~ q) ~\ ' (B - 67) 

(S-'TS^U = — ^ {-m + , r=p + q, (B.68) 



p,q 

where 



x fiip,q) = ffj^n - *s , r=p + q. (B.69) 



For the connected 4-point coupling in leading order we obtain 
(e x ■ e y e x i • ey) c = — — — ^ ^ 



N V 3 ^ i D{p)D{r) D(p')D(r>) 

x {A(q)+X IM ( P ,q)X u {p',q)A flv (q)^ , r = p + q,r' = p' - q , (B.70) 

where 

K ® = wr (R71) 



A^q) H pu (q) - pyH > / + -K 2 b pv = V • (B-72) 



H p (q)H u (q) 1 2 
H{q) 2 pi 

We consider the isotropic case T = L, wo = w\ = w, then in leading order for the 2- 
and 4-point functions: 
V 



wm, 



+ 0{1/N) , (B.73) 



o 



^ ~ ~ - v 2 ^ \ K / \ HJ0)H v (0)~ ' 

x,y,x',y' \ u/ v w 

+ C(l/iV 2 ) , (B.74) 
so that the renormalized coupling, eq. (5.44) is given by 

^ R = ^ {^ (0) + ^o^ (0) ^(0)} + 0(1/*) . (B.75) 
Note (setting to this order mo = M) 

ff(0) = -Ar#o , (B.76) 

*M = , (B.77) 

H ^ = Jv^wrl^ (R78) 

= {^2 V + F 3 (1 - dV)} > (B.79) 
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where 



Hr 



Ly I 

V ^ {p 2 + M 2 ) 2 ' 

1 V f 

V ^ {p 2 + M 2 ) 2 ' 



(p 



2\2 



V V (p 2 + M 2 ) 2 ' 

(p 2 ) 2 - p 4 



(d - 1)V S (p2 + M 2y 



where 



Then 



A^(0) 



w 



\n 2 dw 2 +H 2 - dH 3 



= 7=ri d5 ^ -!} + 



tir 



and so 



A(0) + 



w 



A^(0) = =- + 



\n 2 dw 2 + H 2 - H\/H Q 



w 2 H\ 



H 



o m 



\n 2 dw 2 + H 2 - HI/ H c 



Noting 



Hi = w - M 2 Hn 



H 2 = l- 2wM 2 + M 4 H { 



we get 



Ng R 



+ 



2 (M 2 H - w) 2 



M 2 H m 2 H [1 + \dK 2 w 2 - w 2 /H c 
B.2.2 Lattice artifacts 



0(1/N) 



(B.80) 
(B.81) 
(B.82) 
(B.83) 

(B.84) 
(B.85) 



(B.86) 



(B.87) 
(B.88) 



(B.89) 



In the rest of this section we restrict attention to d = 2 and consider large physical volumes. 
To see the structure of artifacts we expand w and Hq: 



1 „ M 2 , 



M 2 H = ^-^-(2-£)... 



(B.90) 
(B.91) 



4vr 32vr 
where C = -ln(M 2 /32). 

Taking the continuum limit we have k 2 w 2 = q 2 + C(l/L 2 ), with q given by eq. (5.41) 
Using eqs. (B.90), (B.91) in eq. (B.89) (with d = 2) we derive eq. (5.47). 
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B.2.3 Numerical evaluation of #r 

In the infinite volume limit the gap equations are 



and 



where 



J(M) = w , (B.92) 



k 2 (w-j\+M 2 = ^, (B.93) 



(2tt)2 lf(p) + M 2 ./_,, 2vr + 4) 

Inserting w = J(M) into eq. (B.92) one can find numerically M = M(f, k). 



One also has 

Inserting these expressions into eq. (B.89) it is easy to calculate Ng^{aM). 
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